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Abstract 

An initial evaluation is made of extensive 
three-component (3-C) particle image velocimetry 
(PIV) measurements within the wake across a rotor 
disk plane. The model is a 40 percent scale BO- 
105 helicopter main rotor in forward flight 
simulation. This study is part of the HART II test 
program conducted in the German-Dutch Wind 
Tunnel (DNW). Included are wake vortex field 
measurements over the advancing and retreating 
sides of the rotor operating at a typical descent 
landing condition important for impulsive blade- 
vortex interaction (BVI) noise. Also included are 
advancing side results for rotor angle variations 
from climb to steep descent. Using detailed PIV 
vector maps of the vortex fields, methods of 
extracting key vortex parameters are examined 
and a new method was developed and evaluated. 
An objective processing method, involving a 
center-of-vorticity criterion and a vorticity “disk" 
integration, was used to determine vortex core 
size, strength, core velocity distribution 
characteristics, and unsteadiness. These 
parameters are mapped over the rotor disk and 
offer unique physical insight for these parameters 
of importance for rotor noise and vibration 
prediction. 


Symbols 

C roto r blade cho rd , 0. 1 2 1 m 

CV center of vorticity in (x,y) plane, m 

DNW German-Dutch Wind tunnel 

DLR German Aerospace Center 

HART HHC Aeroacoustic Rotor Test 
k grid point indices 

j indices for instantaneous image 

n number associated with analytical 

velocity profile, Eq. 5 
R rotor radius, 2 m 

r c vortex core radius, m 

r ( ' core radius associated with Rankine 

vortex, m 

/ c u c shape factor for vortex velocity profile 

h,v,w velocity components for x,y, z 

coordinates, m/s 

V T rotor hover tip speed, RQ. (218 m/s) 

v c ‘spin’ velocity at rc, m/s 

(x,y,z) TUN windtunnel coordinate system: (x TUN 
positive downstream, y TUN positive 
starboard, z™ positive up). 
x,y,z PIV image frame coordinates, m, Fig. 4. 
a rotor shaft angle with respect to x TUN 

axis, deg 

8 vortex wander parameter, m 

r c circulation within r c , m/s 2 
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a standard deviation of parameters 

derived from instantaneous images 
£1 rotor rotation frequency, rad/s 

co z vorticity normal to 0,;y) plane, s' 1 
xy blade azimuth angle (0° aft), deg. 

Introduction 

Helicopter rotor noise has been measured 
and studied intensely over that last several 
decades in order to understand and predict the 
noise generation as well as determine methods for 
control and reduction. In particular, rotor blade 
vortex interaction noise (BVI) has received much 
of the attention due to its significance during 
descent. When BVI noise occurs, it can dominate 
the radiated noise field and adversely affect 
community acceptance of rotorcraft (Refs. 1-4). In 
addition to BVI noise, rotor broadband noise has 
also been shown to be significant and even 
dominate the radiated noise field for certain 
operating conditions, such as shallow climb (Refs. 
5, 6). The rotor broadband noise sources are 
namely due to blade wake interactions (BWI 
noise), which is the interaction between a blade 
and the turbulent wakes of preceding blades, and 
the interaction between the airfoil blade and the 
turbulence produced in its own boundary layer 
and near wake (self noise). Both BVI and BWI 
noise sources result from the blade interacting 
with the rotor wake and its associated turbulent 
field. In addition to increasing noise, interaction of 
the wake with the rotor also adversely affects 
vibration levels. 

Significant advances in understanding and 
prediction of these noise sources were made in 
part due to the benchmark database made 
available through the HART I test program of 1 994 
(Refs. 7-9). High spatial and temporal resolution 
acoustics, unsteady blade pressures, and blade 
dynamics measurements were obtained for a 40% 
scale model Bo 105 main rotor. In addition wake 
measurements using LDV and LLS were obtained 
only at limited locations (Refs. 10, 11). The main 
objectives were to develop and improve the 
physical understanding and prediction capabilities 
of blade vortex interaction noise with and without 
higher harmonic pitch control (HHC). Significant 
advances in aeroacoustic analysis and their 
validations were achieved, substantially due to 
improved free-wake and prescribed wake models 
(Refs. 12-17). It was determined that HHC 


reduced noise and vibration by increasing blade- 
vortex-miss-distance and altering the wake 
geometry and strength (Refs. 9, 14, 18, 19). More 
detailed wake measurements were found to be 
needed in order to fully understand the effects of 
HHC and to validate the prediction models. 

In the last decade, advances in measurement 
techniques and digital cameras have allowed 
measurements of the rotor wake to be more easily 
attainable. Measurement systems such as 3-D 
LDV, and PIV (Particle Image Velocimetry), are 
now routinely used to acquire instantaneous 
velocity measurements over relatively large areas 
and volumes. Leishman (Refs. 20, 21), Han (Ref. 
22), Coyne (Ref. 23), Bhagwat (Refs. 24, 25), and 
Martin (Ref. 26) have measured the tip vortex 
generated by a one- and two-bladed rotor in hover 
using 3-D LDV. The swirl and axial velocity profiles 
of the tip vortex as it ages as well as derived 
parameters such as core size, strength, and 
turbulence quantities have been documented in 
detail. Martin (Ref. 27) also examined and 
compared the wake (up to one rotor revolution) of 
a single bladed hovering rotor with 4 different tip 
shapes. The initial vortex structure and strength 
was shown to vary significantly as a function of tip 
shape. Heineck (Ref. 28) used 3-C PIV to 
examine the wake of 2-bladed rotor in hover for 
wake ages up to 270 degrees. The effect of 
vortex wander and the effect it may have on 
determining core size from the measured data was 
presented. McAlister (Ref. 29) also used 3-C PIV 
to examine the wake of a hovering rotor. The 
effect of a turbulence-generating device attached 
to the tip of the blade was shown to significantly 
alter the size and strength of the tip vortex. 

PIV measurements have also been made for 
rotors in forward flight. Raffel (Ref. 30) compared 
3-D LDV with 2-C PIV measurements taken at one 
rotor azimuth location. The advantages and 
differences of the measurement techniques for 
application to measuring rotor vortex properties 
were discussed. More extensive 2-C wake 
measurements have been obtained on the 
advancing side of a large model rotor for a number 
or streamwise locations (Ref. 31). Similar 
measurements were also obtained for a proprotor 
(Ref. 32). In that test, multiple vortices with 
opposite sense in rotation were found. 

In 2001, a major international cooperative 
research program Higher-harmonic-control 
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Aeroacoustics Rotor Test II, called HART II (Ref. 
33-36) was conducted under US / German and US 
/ French Memoranda of Understanding, by 
researchers from the German DLR, the French 
ONERA, the Netherlands DNW, NASA Langley, 
and the US Army Aeroflightdynamics Directorate 
(AFDD). HART I is a comprehensive experimental 
program conducted with a 40% Mach scaled 
model of a BO- 105 main rotor in the Duits- 
Nederlandse Windtunnel (DNW) in the 
Netherlands. The objective of the test was to 
obtain extensive 3-C PIV measurements of the 
rotor wake over both the advancing and retreating 
sides of the rotor disk, along with acoustic 
directivity measurements, unsteady blade surface 
pressure measurements, and blade deformation 
measurements in order to (1) to investigate the 
physics of BVI, broadband noise and vibration 
reduction concepts with the HHC technology, (2) 
to develop an analytical prediction capability for 
rotor BVI and broadband noise, and (3) to 
generate a comprehensive experimental database 
for code validation. 

In this paper, an initial examination of the 
HART I 3-C (3-component) PIV wake 
measurements are made in order to determine for 
non-HHC rotor conditions, wake vortex definitions 
for a range of rotor angles from climb to descent. 
PIV vector maps were obtained over the 
advancing and retreating side of the rotor. Details 
of the vortices, such as core size, strength, core 
velocity characteristics, and unsteadiness are 
quantified and determined as a function of rotor 
angle and nominal wake age. Determination of 
these quantities from the measurements using 
different methods are performed and assessed. 
An objective processing method is developed 
and utilized to determine vortex definition. 

Rotor Model Test Description 

The HART II program was conducted in the 
open-jet, anechoic test section of the Large Low- 
speed Facility (LLF) of the DNW, which has an exit 
nozzle of 8 m by 6 m that provides a 19m-long free 
jet with a low-turbulence potential core. 

The set-up for the PIV measurement portion 
of the test is shown in Fig. 1. The rotor hub was 
located nominally 7 m downstream from the nozzle 
exit and maintained at 0.92 m above the tunnel 
centerline. The rotor is a 40-percent, dynamically 


and Mach-scaled Bo105 main rotor with a diameter 
of 4 m. It consists of four hingeless blades that 
have a pre-cone of 2.5 degrees at the hub. The 
rectangular planform blades have a chord length 
of 0.121m, NACA 23012 sections (with tabbed 
trailing edge), -8 degrees linear twist and standard 
rectangular tips. For the HART II test, the rotor was 
operated at a nominal rpm = 1041 (hover tip speed 
of 218 m/s), C T =. 0044, and an advance ratio of 
0.15, for 6 rotor shaft angles, a = -6.9° (climb), - 
3.6°, -0.5°, 2.4°, 5.3°, and 1 1 .5° (steep descent). 
The baseline condition at rotor shaft angle of 5.3°, 
was run for selected HHC conditions. More 
detailed information on the rotor characteristics are 
given in Ref. 36. For this paper, only the non-HHC 
conditions are considered. 

Three-Component PIV Measurements 

The rotor wake was measured on both the 
advancing and retreating sides of the rotor using 
3-C PIV. The measurement locations (cut planes) 
are shown in Fig. 2. The two rotor azimuthal 
orientations, (a) and (b), are used in order to keep 
the blades from interfering with the wake 
measurements. There were approximately 50 
locations on each of the advancing and retreating 
side for the baseline and 2 HHC conditions. For 
limited locations on the advancing side (labeled 
with red numbers) PIV measurements were made 
for the 6 shaft angles. At every the PIV 
measurement location, at least 100 instantaneous 
vector map images were obtained. 

The PIV setup implemented for this test 
consisted of five digital cameras and three double 
pulse Nd:YAG lasers. The lasers and cameras 
were mounted onto a common traversing system 
in order to keep the distance between the 
cameras and the light sheet constant while moving 
to different measurement locations within the rotor 
wake. Figure 1 shows the rotor and traverse 
equipment in the DNW test section. The 5 
cameras were located on the 15m vertical tower 
and the lasers were located underneath the rotor. 
Four of the cameras were for the 3-C PIV 
measurements. The fifth camera was used to 
visually check seeding prior to PIV data 
acquisition. To obtain measurements on the 
retreating side the entire support structure and 
tower was repositioned to the retreating side. 
Measurements could be continued without 
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Figure 1. HART II rotor and 3-D PIV setup in the 
DNW test section. 


recalibration. Further details of the PIV systems 
see Ref. 34. 

At each of the cut planes shown in Fig. 2, 3-C 
PIV measurements were obtained simultaneously 
from two systems. This was done in order to obtain 
both a large image frame of the vortex and its 
surrounding flow as well as a small, higher 
resolution, image frame focused primarily on the 
vortex core region. Figure 3 presents an example 
result from each system (to be subsequently 
discussed). Each camera had a resolution of 1 024 
by 1280 pixels and was digitized to 12 bit. The 
cameras were separated vertically on the tower by 
5.2 m, with one camera from each system above 
and below the rotor plane. For the higher spatial 
resolution system, 300mm lens were used and for 
the lower resolution, larger image frame 
system, 100mm lens were used. 

Three double pulsed Nd: YAG laser systems 
(each 2x320 mJ) were located on the base of the 
common support system as shown in Fig. 2. A 


large (nominallyl.Sm) light sheet of nominally 7mm 
thick was oriented 30.6° with respect to the wind 
tunnel x™ axis. The black cut lines shown in Fig. 
2, indicate the top view of the light sheet 
orientation. The horizontal distance between the 
light sheet and cameras was 5.6 m. The cameras 
and lasers were synchronized with the one-per- 
rev signal given by the rotor, which allowed for 
recording at desired phase-angles of the rotor 
blade. 

Flow seeding was accomplished by a specially 
designed seeding rake located in the settling 
chamber, immediately upstream of the turbulence 
screen. The rake was nominally 3 m x 4 m and was 
connected to Laskin nozzle particle generators. 
Di- 2 -Ethylhexyl-Sebacat (DEHS) was used as the 
seed material. The mean diameter of the particles 
generated was below 1 pm. The seeding rake was 
remotely moved both vertically and horizontally 
during testing to seed the area of interest. 

The procedure was designed to measure the 
wake from the same blade at all cut line locations 
indicated in Fig. 2. This would allow for the tip 
vortex to be investigated in detail from its creation 
through its evolution as it convected and 
traversed downstream throughout the rotor disk 
region. Because of this, any questions of blade- 
to-blade differences could be avoided. (Blade to 
blade differences were however examined in this 
test, but results are not reported here.) 

The measurement process commenced at 
one of the most upstream locations indicated in 
Fig. 2 and progressed downstream along a line of 
constant y-a* (indicated as the black horizontal 
lines). For example, the starting locations for 
locations shown in Fig. 2 (a) are cut lines #1, #9, 
#17, and #24. The reference rotor position is 
identified as vy=70°, which is the azimuth position 
of the blade located in the first rotor quadrant. The 
camera and lasers are synchronized with the one- 
per-rev signal of the rotor and hence PIV image 
data is acquired when the appropriate blade is 
positioned in the reference location. Image pair 
sampling was triggered every 3 rd rotor revolution, if 
the rotor maintained a constant rpm, otherwise the 
trigger could be slightly longer (Ref. 34). Once 
100 image pairs were obtained the common 
support system was traversed downstream to the 
next cut line. The rotor reference blade angle was 
then incremented by 90° in order to continue 
measuring the wake from the same blade. The 
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procedure was repeated in a similar manner for all 
locations shown in Figs. 2 (a) and 2 (b). 


resolution of 0.0058 m / vector. The DLR PIV 
measurement area was smaller (nominally 15.2 cm 
by 12.9 cm) and was centered within the larger 
DNW area on the vortex core region. The DLR PIV 
image shown in Fig. 3 was processed with 
interrogation window sizes of 24x24 pixels, which 
results in a resolution of 0.00154 m/vector, nearly 
4 times that of the DNW data. For the 24x24 pixel 
interrogation window, the maximum vorticity was 
found to be about 16% higher than that obtained 
using data processed with the 32x32 window. 
Velocity vector fields were obtained from each PIV 
image-pair using the cross-correlation method 
(Refs. 34, 37, 38). In Fig. 3 and in presentations 
to follow, the magnitude of velocity vectors for the 
DNW frames, 14.2 m/s equates to a “scale length” 
of 10 cm, and for the DLR frames, 33.3 m/s 
equates to 10 cm. 


a) rotor position for PIV measurements location 
indicated (blade in I s ' quadrant at \|/ = 70°). 


’DNW PIV image frame 


0.24 - 0 . T 6 - 0.08 0.00 0.08 0 . 
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b) rotor position for PIV measurements location 
indicated (blade in 1“ quadrant at vy = 20°). 

Figure 2. Schematic showing the (cut plane) 
locations of the 3-C PIV rotor wake measurements. 


'DLR' PIV image frame 


Figure 3 presents example DNW and DLR 
results for the measurement cut plane location of 
#21 (x™ =0.725 m, y™ = 1.4 m) (see Fig. 2 (b)). 
The DNW obtained PIV measurements over a 
large area that was nominally 45.3 cm by 36.7 cm. 
The PIV images were processed with a 32x32 
interrogation window size that resulted in a 


Figure 3. DNW and DLR PIV measurement image 
frame areas. Every 3 rd vector shown in DNW image, 
every 5 th vector shown in DLR image. 
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PIV data validation 

PIV velocity vector fields can contain incorrect 
vectors (outliers). Before using the PIV data it is 
important to identify and replace the outliers with 
valid data to insure the accuracy of subsequent 
results. Sources of outliers in the HART H data 
included seeding problems, low laser light 
intensity, blade in the field of view (causing 
reflection and obstruction), and large laser pulse 
separation times (Refs. 34). Prior to constructing 
the final 3-C vector fields a scan of the vectors in 
each of the 2-C vector fields are checked. Each 
vector is compared with each of its eight 
neighbors. If the difference in magnitude 
exceeds a specified threshold, its tag is increased 
by one. If a vectors tag value is greater than 5 is 
removed and replaced using a bi-linear 

interpolation from the surrounding valid 

neighbors. The threshold used was always 
slightly conservative in order to ensure that no 
good vectors were removed. 

Vortex Alignment 

For this study, the PIV plane cut 30.7° 
orientation is used for nominal alignment to the 
normal to the vortex axis. Based on pretest wake 
predictions, tilt angles of the vortex with respect to 
the cut lines along y TUN = 1.94, 1.7, 1.4, 1.1, 0.8, - 
0.8, -1.1, -1.4, -1.7, and -1.94 m would be 
expected to be on the order of 0 = 35°, 30°, 20°, 
0°, -3°, 20°, 10°, -3°, -20°, and -35°. Because of tip 
vortex roll-up, there is some variation from this 
along the y TUN cut lines,. In addition, general wake 
unsteadiness can cause instantaneous PIV image 
frame orientation differences. 

It is beyond the scope of this paper to 
determine tilt and account for it in the derivation of 
vortex parameters. A rough estimate is that local 
velocities and dimensions can be in error on the 
order of ±(1-Cos(0)) or about ±15% for 0 = ±30°. 
This is further discussed in the text, but the 
methods used to determine key vortex 
parameters are expected to result in error of less 
than half of this. 

Wake Presentation 

A series of 3C PIV results are presented to 
show the basic character of the rotor wake, as well 


as to successively show derived parametric 
presentations of the wake. 

Figure 4 illustrates the view orientation of the 
PIV cut planes (of Fig. 2 and Table 1), presented 
in the next figures. The view, as well as each cut 
plane is orientated 30.7° from the tunnel axis. The 
DLR and DNW image frames, which coincide with 
the cut planes, are defined with coordinates (x, y), 
where y is aligned with the z™ axis. The frames 
are in rows such that the y axis cuts through the 
same y^ axis value. The (x, y) planes are 
positioned progressively lower (more negative 
Ztun) for increasing Xtu N . This is because the 
camera positions were lowered for the more 
downstream locations, in order to capture the 
vortices. The x axis is offset in angle (30.7°) from 
the downstream (x^ axis) direction. For the 
advancing side cuts, the PIV cameras were on the 
positive yjuN side, out of the flow, so the originally 
processed images were reversed for this 
presentation. For the retreating side, the cameras 
were on the negative y™ side, so no image 
reversal was necessary. 



Figure 4. Sketch to illustrate the view orientation of 
PIV cut planes presented in figures to follow. 
Advancing side cuts are shown. 

Wake Images for Baseline Rotor Condition 

Advancing Side. For a rotor baseline descent 
condition important for BVI noise, Fig. 5 shows 
wake images for six PIV cut planes along y TU N = 
1.4m, viewed in the manner described for Fig. 4. 
Corresponding to Fig. 2 (and Table 1), the cut 
planes are #17, 18, 19, 21, 22, and 23. Because 
the spacing between frames corresponds to Ax TUN 
Cos 30.7° in this view, frames overlap in some 
regions. Where this happens, the frames at larger 
x^n are partially covered. The velocity and vorticity 
data shown in all frames are from PIV images and 
processing by DNW. The outline of the smaller 
DLR frame is shown within each. It contains the 
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primary vortex and is used in subsequent 
determination of detailed vortex parameters. 
Because of its finer spatial resolution, DLR data 
were used to generate all the instantaneous 
centers of vorticity (subsequently defined). 

The top frames of Fig. 5 present processed 
velocity vector maps for instantaneous points in 
time for the wake regions captured. Each arrow 
represent the x and y components of velocity, 
after a mean velocity is removed to permit a (much) 
better viewing of the vortex field. That is, 

U(x,y,z = 0) = ui + vj + wk U = U To , -U 0 , where 
U 0 = uJ+v 0 j+w 0 k\s the mean velocity within the 
region of the larger frame but outside of the 
smaller frame. These U 0 values are listed, in Table 
1 of the Appendix, for the frames shown. (It is 
noted that in previous LDV flow study 8 
presentations for this^ rotor, a set estimated value 
of the equivalent oft/ 0 were used for all locations 
over the rotor disk. But here by this method, one 
sees that there is a substantial variation across the 
rotor disk region - with stronger downwash 
components in the downstream direction.) As 
previously mentioned for the DNW frame 
presentations, the velocity vector magnitudes can 
be equated to the frame size (width is 45.3 cm with 
14.2 m/s equating to about 10 cm). The 
magnitude of the w velocity component is 
indicated by color contour under the ( u , v) 
vectors. Negative w values indicate component is 
directed into the page. For clarity, only every third 
vector along aline is plotted, so only 1/9 m of the 
vectors are plotted. 

The second row of frames from the top of Fig. 
5 contain averaged velocity vectors determined by 
simply summing the 100 instantaneous vector 
maps and dividing by 100. Because of 

unsteadiness (or wander) in vortex position and 
flow details between the instantaneous vector 
fields, smearing of detail occurs, but the basic 
mean wake-shear-flow character is better 
revealed. 

The next two series of frames are of 
instantaneous and averaged vorticity distributions 
generated by the corresponding velocity 
distributions above, that is 

CO z = (Am / Ay — Av / Ax) / 2 . Here Ax and Ay 
correspond to the DNW grid spacing. By the right- 
hand-rule, negative co z indicate rotation is 


clockwise in this view. The vorticity distribution 
presentations, unlike the velocity presentations, 
are independent of the choice of any mean 
velocity values. They are seen to reveal, more 
clearly, the early wake vortex sheets in the 
upstream direction (frames on the left), the roll-up 
evolution, and separation of the blade vortices in 
the downstream direction. For example, one can 
see signs of a near, but not direct, blade vortex 
interaction (BVI) on image frame #22. The lines of 
counter rotating shear from a blade's two 
boundary layers identifies the wake of a recent 
blade passage. This shear-line pair has little net 
circulation and represents the “2D” mid-blade 
wake that eventually dissipates and does not roll- 
up. (A PIV cut plane, if placed nearly normal to 
these present cut planes, would reveal the same 
wake to have net circulation shear lines and thus 
have a tendency to roll-up in that plane (much like 
the view of frame #17). The present data, of 
course, only resolves the PIV cut-plane normal 
component of vorticity.) 

In the bottom two frames of Fig. 5, the centers 
of vorticity, CV, of each of the 100 images (that 
meet criteria listed) are presented. The CV 
positions are used to identify the vortex center in 
each frame for vortex parameter calculations. The 
scatter of these CV positions can also be viewed 
as representing the degree of wander (non- 
repeatability) in the instantaneous images. The 
center of vorticity is defined as 

CV = X(ffl^ •(**,?*)/£ (*> z )* W 

k k 

where k are grid point indices taken over the 
whole DLR frame. For most cut planes, not all 1 00 
instantaneous images have CV positions 
determined. Excluded are images where the 
maximum values of co z do not exceed a 20 % 
threshold of the maximum co z , found by 
considering all 100 image frames. This eliminates 
extraneous images where seeding, laser light, 
other PIV conditions, or an occasional vortex- 
formation disruption produced ill-defined vortices. 
In frame #19 of Fig. 5, for example, 11 of the 100 
images are not included. However, most frames 
only excluded 5 or less. An additional criterion is 
used to properly locate the vortex within each 
frame. In the evaluation of Eq. (1), only those 
values of co, that are of a correct sense (as 
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(a) forward positions (x^ negative) 

Figure 5. Advancing side, Rotor angle 0=5.3° (baseline condition) - velocity and denved parameter maps. PIV image 
frames positioned along y TUN = 1.4 m. 
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(b) aft positions (x TUN positive) 
Figure 5. (Continue) 
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(a) forward positions (X TUN negative) 

Figure 6. Retreating side, Rotor angle a=5.3° (baseline condition) - velocity and derived parameter maps. PIV image 
frames positioned along y T n N = - 1.4m. 
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(b) aft positions (x TUN positive) 
Figure 6. (Continue) 
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determined by evaluation of Eq. (1) for an 
averaged vorticity distribution) and that exceed a 
certain percent of the maximum value of fw z within 
each image are included. In Fig. 5, the results of 
applying 20 % and 80 % criteria are shown. The 
use of 20 % eliminates “noise” within the frame, 
while capturing the vortex and its shear layer. 
Applying an 80 % criterion is found to substantially 
isolate the vortex itself, even when attached to a 
strong shear sheet - such as that seen in frame 
#17 for the young vortex. Because of vortex 
wander the 80 % criterion is seen to produce 
much larger scatter of CV. For some cut plane 
locations, especially where strong shear is present 
and localized “hot spots” occur on instantaneous 
frames, the criteria (especially 80%) can create 
false CV locations in some cases. Although this 
has not been fully evaluated, this is more of a 
problem when the vortex itself is very weak or 
highly fragmented, as would be case for vortices 
after direct BVI. 

Figure 5 gives some illustration of vortex axis 
tilt with respect to the cut plane. Considering 
frame #22, there is not a pronounced w velocity 
component acting through its center, although an 
asymmetry is noted in w across this clockwise- 
rotating vortex in the image frame y direction. This 
indicates that the axis of the vortex (if modeled as 
a rotating tube with axis along the center of the 
tube) is somewhat tilted with respect to the cut 
plane - roughly in the negative x direction from 
the positive z axis (see Fig. 4). Other evidence of 
this tilt is seen in the slight vortex core elongation, 
in contrast to a purely circular pattern - particularly 
seen in the vorticity plots. This elongation is also 
seen in the distribution of CV. 

Retreating Side. Figure 6 has the same 
presentation as Fig. 5 (except for amplitude scaling 
of vorticity, to be discussed) for y^ =-1.4 m on 
the retreating side. The five frames shown are 
numbers #43-47, as identified in Fig. 2. As with 
the rotor advancing side view of Fig. 4, the 
observer is on the retreating side and off the rotor, 
but the view is oriented upstream at 30.7°, rather 
than downstream. As expected, co z is again 
directed into the page, with corresponding 
rotation being clockwise. This indicates that the 
blades that produce these vortices maintain 
positive lift over their tip region, thereby shedding 
vortex / vortex sheets of the same orientation as 
that for y TUN =1 .4 m of Fig. 5. A difference in wake 


character seen for this retreating side is that the 
vortices are stronger and better defined. The 
contour color scaling for vorticity amplitude has 
been extended by a factor of 5 (this can be viewed 
as a dynamic range change) for frames #43, 44, 
and 45, in order to properly visualize the vortex 
core character. (Unfortunately, this tends to wash 
out flow field details.) For frames #46 and 47, a 
factor of 2 was used in order not to lose detail in 
frame #47, but note that this scale somewhat 
unrealistically elongates the core size appearance 
in frame #46. Even at an early age (frame #43), 
the vortex is more fully rolled-up compared to 
advancing side vortices. This would be expected, 
due to higher blade tip loading that concentrates 
the vortex and reduces the extent of the vortex 
sheet to be rolled-up. The vortex wander is 
substantially lower and, unlike the advancing side 
the wander for 80 % criterion, is about the same as 
that of the 20 % criterion. The only exception 
shown is in frame #47 where the vortex is 
fragmented, likely split due to a blade interaction. 
Here, the blade boundary layer trail (pair of 
counter-rotating sheets) is seen to intersect a 
weakened vortex. 

Compared to Fig. 5, a stronger w component 
is apparent over the vortex core region in Fig. 6. 
The ( u , v) components (shown as vectors) are 
correspondingly stronger for these retreating side 
vortices. This w component represents the vortex 
axial component, to the extent that the PIV cut 
plane is perpendicular to the vortex axis. A feature 
similar to Fig. 5 is the asymmetry, across the vortex 
region, in the sign of w. However, the sign is 
reversed (top to bottom) indicating, as anticipated, 
that the vortex is tilted from the cut plane in the 
opposite direction. 

Rotor Shaft Angle Variation - Advancing Side, 

For the advancing side at ynjN=1-4, results 
from PIV cut planes #20, 21, and 22 are 
presented for different rotor shaft angles. The 
angles a = - 6.9° (climb condition), -3.6°, -0.5°, 
2.4°, and 11.5° (steep descent condition) are 
represented in Figs. 7, 8, 9, 10, and 11, 
respectively. These use the same presentation as 
Fig. 5 for a = 5.3°. 

At the largest negative angle, a = - 6.9° (Fig. 
7), the vortices and their sheets are below the 
blade passage region. Here, for the x™ range 
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Figure 7. Rotor angle a = - 6.9°, Advancing side, aft positions (x^ positive) - velocity and derived parameter 
maps. PIV image frames positioned along y TUN = 1 .4 m. 
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Figure 8. Rotor angle a = - 3.6°, otherwise same as Fig. 7. 
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Figure 9. Rotor angle a = - 0.5°, otherwise same as Fig. 7. 












Figure 10. Rotor angle a = 2.4°, otherwise same as Fig. 7 
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Figure 11. Rotor angle a =11.5°, otherwise same as Fig. 7. 
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shown, the blades do not interact strongly with the 
vortices. As the rotor angle increases, the vortices 
and their sheets are seen to approach, 
convergence, and cross through the region 
where the blades interact most strongly. At a = 
11.5°, in Fig. 11, one observes that the blade 
vortices are well above their vortex sheets. At 
intermediate angles a = 2.4°, -0.5°, and -3.6° 
(Figs. 8, 9, and 10), one can see the wake of the 
most recent blade passage - identified, as before, 
by the lines of counter-rotating shear. In frame #22 
for Fig. 10, there appears to be a highly disrupted 
(burst) vortex from a BVI occurrence - much like 
that of frame #47 of Fig. 6. In both cases, CV 
positions are still locatable, although there is 
question as to the presence of false positions (for 
the 80 % criterion CV) or even if one should 
interpret a vortex "per se” as being present. 


Vortex Definition and Discussion 


Processing / Calculation Methods 

Velocity Cuts. Figure 12 shows four velocity 
vector distributions corresponding to PIV cut 
plane #22 for the baseline rotor condition 
presented in Fig. 5. Different instantaneous 
velocity vector plots are shown In Figs. 12 (a) and 
(b) to illustrate a typical degree of variation that 
occurs between different PIV acquisitions. As 
done in Fig. 5, Fig. 12 (c) presents the result of 
straight averaging of all 1 00 instantaneous fields, 
which renders a smoothed velocity field 

presentation. As previously mentioned, this type 
of averaging can result in smearing of the vortex 
flow details - because of vortex wandering. In 
each of the instantaneous plots, Figs. 12 (a) and 
(b), the respective CV locations (determined using 
the 80 % criteria) are shown. In Fig. 12 (c), all 1 00 
CV locations are shown to indicate the spread or 
degree of wander as presented in Fig. 5. Figure 
12 (d) shows the result of conditional averaging, 
which is accomplished by first aligning CV 
locations prior to averaging. This is seen to 
improve the sharpness of the vortex definition. 
Note that the image frame is reduced in size due 
to the removal of non-overlapping regions 
occurring when the CV locations are shifted. 

Figure 13 shows vertical and horizontal 
velocity cuts taken through the CV for the 


instantaneous and averaged distributions of Fig. 
12. For Fig. 13 (c), an average CV position is 
used. The horizontal velocity components, 
plotted along the vertical cut, and the vertical 
components, plotted along a horizontal cut, give a 
swirl appearance. Upon defining v c to equal 1/2 
of the difference in peak velocities (identified in 
figures) across the vortex, and r c to equal 1/2 of 
the distance between the locations of these peak 
velocities, the horizontal and vertical cuts give (r c , 
v c ) values of (.0139m, 5.20m/s) and (,0169m, 
6.54m/s), respectively, for Fig. 13(a). For Fig. 
13(b), these are (.0108, 5.19) and (.0299, 6.05) 
with units of m and m/s. The ragged character of 
the instantaneous image cuts give uncertainty to 
“true accuracy” of the choice of r c andv c . The 
magnitude of v c is always biased to higher values. 
Figure 13(c) renders ( r c , v c ) values of (.0299, 
3.37) and (.0475, 4.23), which does illustrate that 
simple averaging produces a larger r c and smaller 
v c compared to instantaneous values. For Fig. 
13(d), the values are (.0185, 4.09) and (.023, 
4.67), which should be closer to accurate values. 
But still, the results are still biased because 
individual vortices of different sizes tend to smear 
and flatten upon averaging - even if centered. 



Figure 12. Velocity vector distributions over image 
frame #22, for a = 5.3° showing instantaneous and 
averaging methods. 
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a) image # 5 


b) image # 30 
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Figure 13. Vertical and horizontal velocity cuts of 
the respective Fig. 12 distributions. 


vortex than that found above using the two cut 
procedure of Fig. 13 (a) and then averaging the 
two results. The analytical curve in Fig. 15 (b) is a 
calculated description of the velocity profile using 
a method to be described. 



u.vw 

x[m] 


Figure 14. Integration disk with radius r to determine 
circulation T. 


Vorticitv Disk and Analytic Velocity Profile 
Method. An alternate method is developed with 
the intent to determine r c and v c that are better 
based statistically, less bias, and in a form that 
renders a more complete description of the 
vortices. Instead of operating on a limited number 
of velocity cut distributions, vorticity distributions 
are used. A present limitation to be mentioned is 
that the procedure intrinsically assumes that the 
vortex tube is aligned normal to the cut plane. 
Errors associated with this are subsequently 
discussed. 

Figure 14 shows the vorticity distribution 
associated with the instantaneous image of Fig. 
12 (a). The CV position is shown with a radius r 
defining a circle (or disk) that encloses vorticity. 
The procedure of the “disk” method is to first 
integrate (both positive and negative) vorticity 
over areas within succeeding values of radius r, to 
render circulation r asa function of r. The result 
of this integration is shown in Fig. 15 (a). Using the 
relationship v = T/2Kr, one obtains Fig. 15 (b), 
where r c is chosen at the maximum value of v 
(designated asv c , v c = r c Hitr c ). Here, this 
method gives (r c , v c ) = (,0143m, 4.9m/s), which is 
deemed to be a more accurate measure of the 


a) circulation 


b) swirl velocity 



Figure 15. Determination of vortex parameters, (a) 
circulation and (b) swirl velocity. 



Figure 16. Types of velocity profile shapes 
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Besides providing values of r c and v c , the 
disk method velocity profile in Fig. 15(b) can be 
presumed to be a good representation of that for 
the vortex. As a first step in quantifying the profile 
shape of the vortex core region, Fig. 16 illustrates 
several possible profile shapes within the core of a 
vortex. Convex, linear, and concave profiles are 
shown. If A in the figure is the area of velocity 
deficit, then A = ( l/2)/ c v c , where / c is the 
Rankine (linear profile) vortex core radius with 
velocity v c . In general, for the area 

A = v c r c \^(\-vlv c )dr (2) 

where r= r/r c , so 

/ c /r c =2j' o (\-v/v c )dr (3) 

The ratio r’ c /r c is taken to be a “shape factor 1 ' 
of interest. When r' c !r c < 1 the profile is convex, 
when / c !r c = 1 the profile is Rankine (linear), and 
when i'q lr c > 1 it is concave. Equation (3) is 
readily evaluated from processed v results, such 
as in Fig. 15 (b). When / c !r c < 1 , we can chose an 
algebraic velocity-profile description to describe 
most of the profiles found in the present study. 
This is the Vatistas 39 profile, where the velocity is 
given by 

v/v c = 2 1>n {? /(l + r 2n ) l,n ) (4) 

where within the vortex core, ? <1 , and outside of 
the core, r>1. The number n= 1 gives the Scully 
vortex model 40 shape, often used for rotorcraft 
vortex definitions 41 ; and n= 2 gives a vortex 
description 22 ' 24 found to give good results for 
vortices measured for a model rotor in hover. The 
Rankine vortex model is obtained when n =«>. 
Equation (4) does not provide a description for a 
convex profile. From Eqs. (3) and (4), one can 
relate values of n (such as .5, 1, 2, and °°) to 
values oi/ c /r c (.455, .614, .754, and 1, 

respectively). 

The analytic curve of Fig. 15(b) is determined 
from Eq. (4), using the r c , v c , and n values found 
in the use of Eq. (3) with the data. The same 
image #5 results are presented in Fig. 17 (a) in the 
swirl format of Fig. 13 (a) for direct comparison with 
the velocity cuts. Also, Fig. 17 (b) compares with 


Fig. 13 (b) in the same way, but for image #30. 
Note that a very different value of n was found for 
this particular velocity profile. 

It is important to note that emphasis is given in 
this paper to define velocity profile characteristics 
within the vortex core. In Fig. 17, the good 
agreement of Eq. (4) using the respective n values 
is not matched outside of the core. The n values in 
this equation determine how much of the 
circulation associated with the vortex is contained 
within the core (for the particular equation form). 
With T c being the circulation within the core (as 
per the disk method, v c = r c /27tr c )and r T0T _ Vortex 
is the total circulation associated with the vortex, 
then 

I~c = 2 1/n r TOT _ VorJa (5) 

where for«= 1, 2, «, 2~ 1/n = .5, .707, 1. It is seen 
that because of Eq. (5), one could emphasize the 
outer velocity profile in determining n, to get better 
agreement there, but again the present purpose 
is to describe core region profile. 

Figure 18 presents the analytic result for the 
same PIV cut plane #22 using average ? c , v c , and 
n values (.0183m, 4.69m/s, and .684, 

respectively). These were obtained from simple 
mean averaging of r c and v c , and for n, 
determined from corresponding mean values of 
/ c ir c . This result can be taken to represent the 
mean vortex of the PIV cut plane. The parameters 
r c , v c , and / c !r c (not n) were found to have near 
normal statistical distributions. This averaging 
method gives less bias error than one involving 
conditional averaging of profiles and then 
determining parameters (see previous section.) 
To address unsteadiness (non-repeatability) 
questions for instantaneous images, Fig. 18 
“unsteadiness” limits are shown that enclose the 
mean vortex velocity profile. The limits provide a 
measure of how well an instantaneous vortex is 
represented by the mean vortex. This is not an 
“error” band. The limits enclose curves 
determined when r c , v c , and n are each 
independently varied about their standard 
deviations, ±a Tc , ±cr Vc , and ±a n (determined from 

/ c lr c ±o r ^i Tc )■ 

As mentioned, the above method assumes 
that the vortex axis (tube) is normal to the cut 
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Appendix 


Table 1 . Wake vortex parameter definitions for 80 % CV criterion 
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Several important determinations and 
confirmations for wake vortices were found. In 
general, previous predictions and measurements 
(by other, less accurate means) confirm the vortex 
strength dependence with blade tip loading. 
However, vortex core size can be substantially 
smaller (-40% of previous values on the 
advancing side near BVI locations). Also, previous 
assumptions about core velocity profile-shape 
definitions (such as Scully’s model) can be 
substantially in error. Large ranges of vortex 
velocity shapes exist over the rotor disk for this 
forward flight case. Unsteadiness or wander of 
vortex path starts out small at creation but grows to 
about 15 to 20 % of the chord length. 
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Figure 21. Core size variation 
due to rotor angle change. 
Advancing side. CV criterion of 
80%. 



Figure 22. Swirl velocity 
variation due to rotor angle 
change. Advancing side. CV 
criterion of 80 %. 
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Figure 23. Vortex wander 
variation due to rotor angle 
change. Advancing side. CV 
criterion of 80 %. 


and 11.) As in Fig. 19, the centers of vorticity were 
chosen using the 80% criterion. The basic tends 
are, as the rotor condition changes from one of 
climb (a = -6.9°) to steep descent (a=1 1 .5°), that 
the vortex size increases on the order of 1 00 %, 
the vortex swirl velocity decreases about 60 %, 
and wander increases about 50 %. 


The 3-C PIV database from the HART I 
program, with detail flow information over the rotor 
disk for a rotor in forward flight, is unique and 
provides a means to evaluate and develop 
improved wake codes. The present effort was the 


first in-depth look at the data. The present set of 
results can give substantial insight to the rotor 
wake’s vortex structure. 

Processing methods to define important 
vortex parameters from the PIV database were 
examined and a new method was developed and 
evaluated. Compared to commonly used vortex 
definition methods using velocity cuts of 
instantaneous and/or averaged (straight or 
conditionally) vortex velocity fields, the present 
center-of-vorticity criterion and vorticity “disk” 
integration method appears more objectively 
based and is robust. (As with any method, 
individual results should still be subjected to 
interpretation query.) 
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(b) CV criterion of 80%. 


Figure 20. (continued) 
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taken as being more correct, based upon 
comparison with Fig. 5 and other data. The 
velocity contours show only a slight decrease with 
CV 2 o from those of CV 80 . The shape factors show, 
in general, an increase due to CV 20 locations 
being somewhat displaced from maximum vorticity 
locations. Here, the disk method would 
emphasize velocity growth at larger r values, 
thereby rendering more concave profile shapes. 
The wander 5 of the CV position is seen to be 
less with CV 20 . This indicates that, within each 
instantaneous image frame, the region of strong 
vorticity is more stable than the vortices 
themselves, whose instantaneous positions are 
better reflected by CV 80 . 

Upon accepting CV 80 as giving more accurate 
parameters, it is seen that younger vortices 
(negative x TUN near the rotor disk edge) start out at 
with radii r c of 5 to 10 % of chord. The vortices, at 
mid y-ruN values, experience moderate radii growth 
to 10 to 15 % of chord and then appear to start 
contracting. In Fig. 5, some tightening of the 
vortices are seen at positive x^m, but also the 
vortex experience splitting and the disk method 
with CV 80 may have picked up on just one of the 
split vortices. The velocities v c undergo a slight 
drop after the vortex is formed and then remains 
stable. 

Strong vortices are found near vy » 90° (y-ruN = 
2.0). As \y approaches 180° (y-ruN =0. x-rjn is 
negative), they weaken. This is anticipated from 
predictions 13,14 and previous measurements for 
this rotor system 9 . The blades have high tip 
loading near vy = 90° but reduce almost linearly to 
vy = 180°. The fact that r c and v c values remain 
almost constant across much of x^ shows that 
once formed, the vortices are stable and generally 
maintain their strength. Although no detail 
comparative evaluation has been made, the 
present more-accurate PIV vortex measurements 
appear to find the vortices to be much smaller, 
about 40% of that found by the previous LDV 
measurements 18 of this rotor. The data sampling 
for the LDV method in many respects duplicated 
the spatial averaging effects illustrated in the 
straight averaging approach of Fig. 12 (c) of the 
present paper. For that illustration, the straight 
averaging approach produced a radius of about 38 
mm, Ref. 9 found about 30 mm for a similar 


location in the wake, and the present disk method 
finds 14 mm. 

Over most of the advancing side, where the 
shape factor i' c lr c < 1 , the vortices have the 
traditional convex velocity profile shapes. For 
larger y™, values, the commonly assumed Scully 
shape (r£ lr c ~ .61, n = 1) and the more peaked 
profiles of n * 2 ( / c !r c ~ .75) are found. However, 
more inboard, very sharp profiles approaching the 
Rankine profile ( /r c =1) are found. 

The term 5 represents the degree that 
vortices wander within the data set. It is simply 
defined as the average distance 8 that a vortex, 
positioned at CVj within an instantaneous image 

j, is from the mean vortex position, CV , for that 
set of 100 images. 

<5=-Y|cv,-cv| (6) 

j j 

It is seen in Fig. 19 (b) that 8/C starts at negative 
x TUN to be small and increases somewhat in the 
downstream direction. 8 is larger (up to 20% 
chord) for weaker or disrupted vortices. 

Retreating Side - Baseline Condition. Figures 
20(a) and (b), for the retreating side, show that in 
most respects the results for both CV m and CV 20 
are very similar, showing the vortices to be 
generally well defined. The very strong vortices in 
evidence are downstream of high blade tip loading 
as cited in Ref. 9 and 14 for this rotor. At younger 
stages, ? c is small at about 4 % of chord length. 
Downstream, it appears to increase by a factor of 
3. Swirl velocity correspondingly drops. This is in 
contrast to the advancing side where the (weaker) 
vortices appear to stabilize. The shape factor 
/ c !r c indicates that the velocity profiles fit the 
Scully or flatter (rounder) peak profile models. The 
wander 8 is seen to be quite small for young 
vortices and progressively increase to about 15 % 
of chord downstream 

Rotnr Angle Variation- A dvancing Side. For 
variation in rotor shaft angle, Figs. 21, 22, and 23 
show the behaviors of r c /C, v c /V T , and 8/C, 
respectively, over a portion of the first rotor 
quadrant. (One can correspond the grid points 
across yruN=1-4 m to image frames in Figs. 7, 8, 9, 
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a) image #5 b) image # 30 
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Figure 17. Velocity cut presentation of Disk method 
results to correspond to Fig. 13 (a) and (b). 



Figure 18. Mean vortex velocity profile for PIV cut 
plane #23, for a = 5.3°. Unsteadiness limits of 
instantaneous values are shown. 

plane. As previously indicated, the vortices of this 
study can have tilt angle with respect to the normal 
on the order of 30° at the larger y TUN value cut lines 
(outer part of rotor disk). This suggests possible 
local errors of ±15% in velocity and vortex 
dimensions. However, because of an averaging 
effect, the disk method for such a tilt amplitude 
should render a result for r c that is about 7% too 
large and a result for v c that is about 7% too small. 
Such errors should not affect the following 
parametric presentations significantly, but this 
should be kept in mind. 

CV Criterion using Disk Method. The disk 
method is an objective means to determine vortex 
parameters. What it calculates depends on the 
choice of CV in the flow field. The use of the 80 
% criterion, CV 80 , may correctly pinpoint the 


center of the dominant vortex within the image, or 
may locate a “false” location when the vortex is ill 
defined, such as when it’s weak, or it may identify a 
portion of a burst or fragmented vortex. In the 
latter case, the calculation may find a smaller vortex 
than it should. Such results for vortex radius, 
velocity, and shape factor would be averaged into 
the other image results to get a mean vortex. The 
use of the 20 % criterion, CV 20 , tends to find the 
group center of vorticity. Here the CV would 
generally be between areas of peak vorticity and 
thus would tend to calculate a larger dispersed 
vortex. This may be appropriate for burst vortices, 
for example, but it may not be for a small vortex 
attached to a strong elongated shear layer. For 
well-defined vortices, the results for the different 
criteria should be very similar. This in itself makes 
comparison of results from the different criteria 
valuable. 

Vortex Parameter Results 

In this section, vortex parameters are 
examined over the wake for the baseline a = 5.3° 
rotor condition and a variations. The normalized 
parameters are core radius r c /C, mean maximum 
spin velocity v c IV T , velocity profile shape factor 
/ c /r c , and vortex wander 5/C (to be defined) 
Here, V T is the rotor tip speed and C is the blade 
chord length. Many of these parameters are also 
listed in Table 1. 

In the following parameter presentations for 
the baseline case, the results of both CV criteria 
are shown. 

Advancing Side - Baseline Condition. Figures 
19 (a) and (b), using the CV 20 and CV m criteria, 
respectively, show parameter contours plotted 
over the rotor’s advancing side. These are top 
views that show the rotor plane disk edge and 
rotor hub axis defined with respect to the tunnel 
coordinates x^ and y TUN . The centers of the PIV 
image frames are shown as grid points for the data. 

The different criteria are seen to give very 
similar results over much of the rotor disk, 
indicating that such vortices are well defined. As 
expected, based on the above discussions, the 
CV 20 results have larger r c over regions where 
young vortices are forming and where the vortices 
are weak or fragmented. The CV 80 results are 
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